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ABSTRACT 

We explore the cosmological halo-to-halo scatter of the distribution of mass within 
dark matter halos utilizing a well-resolved statistical sample of clusters from the cos- 
mological Millennium simulation. We find that at any radius, the spherically-averaged 
dark matter density of a halo (corresponding to the "smooth-component") and its 
logarithmic slope are well-described by a Gaussian probability distribution. At small 
radii (within the scale radius), the density distribution is fully determined by the 
measured Gaussian distribution in halo concentrations. The variance in the radial dis- 
tribution of mass in dark matter halos is important for the interpretation of direct and 
indirect dark matter detection efforts. The scatter in mass profiles imparts approxi- 
mately a 25 percent cosmological uncertainty in the dark matter density at the Solar 
neighborhood and a factor of ~3 uncertainty in the expected Galactic dark matter 
annihilation flux. The aggregate effect of halo-to-halo profile scatter leads to a small 
(few percent) enhancement in dark matter annihilation background if the Gaussian 
concentration distribution holds for all halo masses versus a 10 percent enhancement 
under the assumption of a log-normal concentration distribution. The Gaussian nature 
of the cluster profile scatter implies that the technique of "stacking" halos to improve 
signal to noise should not suffer from bias. 
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1 INTRODUCTION 

Overwhelming evidence indicates that most of the mass in 
the Universe is composed of dark matter. In the concor- 
dance A Cold Dark Matter (ACDM) cosmological model, 
structure formation is dominated by the gravitational evo- 
lution of dark matter. At the present epoch, most of the mass 
has been assembled into self-bound halos, which are hosts 
to the galaxies, clusters, and groups that are observed. Al- 
though the exact nature of dark matter is unknown, theoret- 
ically motivated extensions to the standard model of particle 
physics suggest cold dark matter candidates which were in 
thermal equilibrium in the early Universe and interact only 
weakly with baryonic matter. 

Weakly Interacting Massive Particles (WIMPs) are 
ideal dark matter candidates as they arise naturally in many 
extensions of the standard model of particle physics. The 
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strength of their interactions can mimic the physical be- 
haviour of the dark matter (weak as well as gravitational) in- 
ferred from a broad range of observations. WIMP dark mat- 
ter candidates h ave a small but non-zero cross-section fo r 
self-annihilation jjungman. Kamionkowski fc Griest lll996l V 
Indirect detection experiments look for the by-products 
of this annihilation, typically in the form of high-energy 
photons, neutrin os and posi trons, a s well as low-energy 
antiprotons (see Jungman. Kamionkowski fc Griest 1 1 19961 : 



iBertone, Hooper fc Silk Il2005ll 

In recent years, halo structure has been widely ex- 
plored using cosmological numerical simulations. Pioneer- 
ing work by Navarro, Frenk & White (1996, 1997; NFW 
hereafter) used numerical simulations to show that the 
spherically-averaged radial density profile of dark mat- 
ter halos is approximate l y "un iversal" . However, NFW 
and later iBullock et al. I (|200ll ) noted significant varia- 
tions in the profile between different halos, in that some 
are better fit by steeper (higher concentration) profile 
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forms than others. These variations wer e shown to be 
corre lated with halo formation time (e.g. IWechsler et al. I 
Upon close inspection, the profile of any individ- 
ual halo cannot be described by any particular smooth 
functional form (e.g. [Reed et al. 1120051: icao fc White II2OO6I : 



iKnollmann. Power fc Knebe I I2OO8I : iLukic et al. I l2009l l. In 

this sense, the halo mass profile is not truly "universal". 
Thus, due to this fundamentally non-smooth nature of the 
radial mass distribution of halos, functional forms do not 
provide a complete and accurate description of the halo den- 
sity. 

In this paper, we quantify the mean and scatter of the 
halo radial mass distribution without the a priori assump- 
tion of a smooth functional form. We use a large sample 
of halos extracted from the "Millennium" cosmological nu- 
merical simulation, whose combination of large volume and 
fine mass resolution are ideal for studying statistical vari- 
ations in halo mass profiles. Our approach begins with an 
empirical non-parametric measure of the distribution func- 
tion of halo densities, allowing us to include the effects of 
halo-to-halo and intra-halo non-universality in addition to 
effects relating to scatter in halo concentrations. We focus 
on small radii where the "smooth" dark matter component 
dominates the mass distribution; here, self-bound satellite 
"subhalos" are deficient bec ause of efficient tidal stripping 
before reaching small radii jSpringel et al. II2OO8I ). This al- 
lows us to quantify the cosmological scatter in the smooth 
component of the dark matter density. We apply our results 
to several cosmological applications, including direct and in- 
direct experimental efforts of dark matter detection. 

Annihilation of dark matter scales as the square of the 
number density of particles, and thus any detected annihi- 
lation signal will be sensitive to the precise distribution of 
dark matter. Past work focused on using extremely high res- 
olution simulations of individual dark matter halos of Galac- 
tic mass to estimate ann ihilation luminosities for particular 
dark matter candida tes (|Kuhlen. Diemand &: Madau 11200^ : 
ISpringel et al. 1 120081 ) . These impressive numerical simula- 
tions were able to quantify the level at which substruc- 
ture contributes to the annihilation signal, as well as get a 
glimpse of the phase-space structure of the Milky Way halo 
at the Solar neighbourhood. However, halo-to-halo varia- 
tions in the radial mass profile implies a "cosmological" un- 
certainty in the predicted annihilation rate in halos, which 
is of course in addition to the uncertainties related to the 
mass, cross-section, and other properties of the dark matter 
particle. 

Dark matter annihilation in halos produces a cos- 
mological background whose strength depends upon the 
numbers of halos throughout the his tory of the Uni- 
verse and their density profiles (e.g. luilio et al71 |2002| : 



IZavala. Springel fc Bovlan-Kolchin 1120101 ) . When integrated 
over all halos, a scatter in halo density profiles implies a 
boost of the annihilation background with respect to the 
case wh ere all halos follow a universal profile without scat- 
ter (e.g. lUllio et"aril2002l ). The strength of the annihilation 
background will thus be sensitive to both the average and 
the scatter in halo dark matter profiles. Inferring the dark 
matter particle mass and cross section from a background 
annihilation signal will thus require separating out the inte- 
grated effects of cosmological scatter. 

In addition, even within the Milky Way, interpretation 



of both direct and indirect dark matter detection efforts re- 
quires an understanding of the intrinsic variance in halo pro- 
files. The cosmological uncertainty in the local dark matter 
density at the solar radius is essential in the interpretation 
of a detection (or lack of) in direct detection experiments. 
In particular, knowledge of the intra-halo scatter in density 
with radius is needed to evaluate results of direct (local) 
dark matter detection efforts in the context of indirect (non- 
local) detection experiments focussed on the Galactic center 
or elsewhere. Thus, both direct and indirect dark matter de- 
tection efforts face the challenge of disentangling the influ- 
ences of cosmological halos from the properties of the dark 
matter particle. 

In this paper, we use a large number of cosmological 
halos to quantify the cosmological variance in halo densi- 
ties, and to assess the level at which the distribution of halo 
densities affects the interpretation of dark matter direct and 
indirect detection experiments. In §[21 we review the general 
characteristics of dark matter halo mass distribution, and we 
give an overview of the cosmological simulations that we use 
in § O We present our analysis of the mass distribution and 
its scatter within halos in the simulation in §|4l showing that 
the distribution of halo concentrations yields an accurate de- 
scription of the measured scatter of the smooth component 
of dark matter within halos. We apply our findings to exper- 
imental searches of dark matter in §(5] we discuss limitations 
of our work (§IU, followed by a brief conclusion (§[?}. 



2 THE DARK MATTER HALO PROFILE 

As a baseline reference for examining the cosmological scat- 
ter of halo density profiles, it is convenient to parameter- 
ize the density profile by a sp herically-averag ed functional 
form. Recent works favor the lEinasto lll965l) profile form 
as a d escription of cosmological halos "('Nav arro et al. II2004I : 
iGao e t al. 2008; Hayashi & White 2008). In this case, the 
logarithmic slope of the density is a simple power-law: 



din p 
dlnr 



= -2 



r^2 



(1) 



The "scale radius", r = r_2, has a density slope of —2, and 
defines the halo concentration as: 



Vvir r20Q 
Cvir — , or C2OO = , 

r-2 r-2 



(2) 



where in the first definition, r^i,- is the virial radius of 
the halo, defined as a sph ere of 95.4 times critical density 
(jEke. Cole, fc Frenk |[l99^ ). while in the second definition, 
r2oo is the radius where the enclosed density is 200 times the 
critical density of the Universe. The parameter a in Eq. [1] 
varies weakly with mass and redshift. On average a = 0.19, 
for z = clusters over th e mass range that we explore (as 
shown in lGao et al. Ilioosl ). Density is given by: 



p{r) 



P-2 exp 



r 

r-2 



(3) 



The normalization of the profile is obtained from requiring 
that the mass of the halo is M = J^^'^" 4:np{r)r^dr, thus 

23/-ciooM 



P-2 = 



a ' a 



(4) 
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Figure 1. Left panel: The probability distribution function (PDF) of concentrations for the 3501 halos with M ^ 8.6 X vS^^h~^ Mq 
shown against best-fit Gaussian with a mean (c) = 4.55, and standard deviation CTc = 1.28. This provides a better match than the best-fit 
log-normal distribution with a mean (logjg c) = 0.669 and a standard deviation in logc of (Tiog^, = 0.121. Concentration fits made to the 
NFW profile produce a modestly narrower distribution, even though the Einasto profile is overall a better description of the halo mass 
distribution. Right panel: The distribution of concentrations in a narrow virial mass range of 8.6 X 10" s£ [M/Zi-IMq] s: 1.7 x 10^"'. 
The shape of the concentration PDF remains Gaussian. 



Here, T{x) and T(x, y) are the Gamma, and Incomplete 
Gamma functions respectively. 



3 SIMULATIONS AND HALO CATALOGUE 

We utilize the Rravity-only N-bod y particle Millennium Sim- 
ulation of Spring el et al. I l|2005l ) . This simulation evolves 
2160^ particles in a periodic box of 500/i~^Mpc usin g 
the gravity solver code lgadget2 (jSpringel et al. Ilioosh . 
a modified ver sion of the publicly available GADGET2 
ISpringel I |2005| ). The cosmological parameters used are 

= 0.25, = 0.75, hubble constant h = 0.73, ^It = 
0.045, n = 1, with power spectrum normalization erg — 0.9. 
The matter power spectrum u sed to create initial conditio ns 
is produced using CMBFAST l|Seliak fc ZaldarriagaHligQe !). 

For this study, it is important to use as large a statis- 
tical sample of halos as possible. In addition, these halos 
must be resolved in the innermost regions in which we are 
interested. We thus consider all halos with more than 10^ 
particles, corresponding to Mvir ^ 8.6 x 10^'^h~^ Mq. This 
results in 3,501 halos. The density profiles for these halos are 
resolved down to ~ 1 — 2%rj,ir, based on c onverg ence tests in 
iMoore et al~l (|l998l '). iPower et al. I (|2003l '). and iReed et~an 
( 20051 ) for halos with similar numbers of particles. 

We define halos an d their cent ers using th e same pro- 
cedure as was done in iNeto et aL ] (12007, ) and iGao et al. I 
(|2008h . Halos are identified initially using friends- of- friends 
with linking length of 0.2 times the mean inter-particle spac- 
ing. Halos are centered on the location of the deepest poten- 
tial of the main subhalo. Halo mass is then determined by a 



sphere of 95.4 times critical density (Mvir), and additionally 
by a sphere of 200 times critical density (M2oo). We deter- 
mine a spherically-averaged logarithmically binned density 
profile for each halo. 



4 COSMOLOGICAL VARIATIONS IN THE 
HALO PROFILE 

In this section we discuss the cosmological variance in the 
properties that describe the profile of dark matter halos. 



4.1 Halo Concentrations 

Halo concentrations have been shown to have significant 
halo-to-halo scatter, with a median that decreases with 
increasing mass a nd redshift (e.g . ^uUock et al. 200 l|; 
iNeto et al. I l2007l : iGao et al. I 1200^ : iMaccio et alT I |2008| ). 
Early work suggested that the distributio n of halo concen- 
trations is log-normal (|Bullock et al. |[200ll ). However, larger 
samples of higher resolution halos reveal significant depar- 
tures from a log-normal scatter, primarily due to a tail of low 
concentrations inferred fr om "unrelaxed" halos (|Neto et al. I 
I2OO7I ; iMaccio et al. I|20C)8|). which tend to conform poorly to 
smooth functional fits ( Lukic et al. |[2009l ). In fact, the dis- 
tribution of halo concentrations is very well descri bed by 
a simple Gaussian, as noted by iLukic et al. I (|2009i ). when 
all (relaxed and unrelaxed) halos are considered. In Fig. [1] 
we show the concentration probability distribution function 
(PDF) from our full sample. We find that a Gaussian de- 
scription of concentrations is a better fit than a log-normal 
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distribution to the halos in our sample. More complicated 
functional d escriptions of halo concen trations such as that 
suggested bv lNeto et aTl l|2007l ') or bv lMaccio et al. I (|2008l ) 
appear unnecessary to describe our data (which consists of 
the high mass halo s ubset of the halos of iNeto et al. ||2007| 
and iGao et al. |[200^ ) . 

In estimating halo concentrations, we consider both the 
Einasto profile (Eq. and the NFW profile (NEW 1996; 
1997): 



(r/r_2)(l + r/r_2)2- ^' 

A concentration defined by the Einasto profile is, in prin- 
ciple, equivalent to that defined by the NFW profile, both 
having a scale radius at r_2. Ifowever, because simulation 
halos tend to better match the Einasto form, which is steeper 
than NFW at the smallest radii, a concentration inferred 
from the NFW form can be biased high or low , depending 
on the range of radii used in fitting (see iGacTet al. 2008. ). 
We find that although the Einasto profile produces a bet- 
ter fit to stacked halos, the distribution of concentrations 
is modestly narrower when fit according to an NFW pro- 
file (see Fig. [1]), with approximately the same mean value 
{c2oo,NFW = 4.55 versus C200, Einasto = 4.50). The Einasto 
distribution remains wider, whether or not the Einasto pa- 
rameter a is fixed or allowed to float as a free parameter 
in the fit. For this reason, we determine halo concentrations 
by fitting the NFW profile (Eq. [5]) in the remainder of this 
paper. Profiles are fit to logarithmic radial bins over a range 
of 0.05 — lr2oo with normalization set by the mass contained 
within r2oo. 

We have confirmed also that fixing the density normal- 
ization instead of allowing it to float as a fit parameter has no 
significant effect upon the derived concentrations. As a fur- 
ther test, we show that with a narrower sample mass range 
of 8.6 X 10^^ ^ [M^ir/h-'^MQ] ^ 1.7 x lO", a Gaussian 
concentration distribution is still preferred over log-normal 
(right panel of Fig. [T]). The similar shape is not surprising 
because the mass dependence on mean halo concentration is 
relatively weak, and more importantly, this shows that the 
concentration distribution is not sensitive to our choice of 
mass range for the halo sample used throughout the paper. 

4.2 Halo Densities 

It is important to note that the Einasto profile has 
been found to fit halos we l l for a "stacked " ensemble 
(e.g. iHavashi" fc White I [200i : iGao et al. I bOOSl ). However, 
the presence of substructure and other peculiarities implies 
that any particular halo profile tends to have significant vari- 
ations from this mean smooth function. 

In the left panel of Fig. [2l we show this scatter in the 
density profile for 100 random halos. Large scatter is also 
seen in the logarithmic slope of the density profile when 
plotted versus r/r_2 (right panel of Fig. [2|. This represen- 
tation removes differences that arise due to concentration. 
The mean slope of the complete halo sample is well described 
by an Einasto profile, and is poorly matched by an NFW 
profile. 

In Fig. O we show the probability distribution function 
(PDF) of densities at various radii in spherically-averaged 
shells from the halo sample. The width of the halo-to- 



halo scatter of density decreases toward larger radii. A 
possible explanation for this radial trend results from the 
fact that the central structure of the halo is assembled 
at higher redshifts than the outer parts of the halo (se e 
e.g. iFukushige. Kawai. fc Makino II2004I : iReed et al. 1120051 1. 
If one assumes that halo density at a particular radius cor- 
relates with the mean density of the universe at the time of 
mass infall, then scatter in mass assembly redshift from halo 
to halo would yield density variations that would be larger 
nearer the center due to the (1 -|- z)'^ evolution of the mean 
matter density. Note that, at all radii, the width of the dis- 
tribution is small compared to the statistical measurement 
uncertainty, which is estimated from Poisson counting of 
particles in radial bins. 

The halo density PDF is well-described by a Gaussian 
at each radius. As an example, in Fig. 3] we show the den- 
sity distribution function at a radius of r = 0.03r2oo. This 
scatter in densities is primarily due to the distribution in 
halo concentrations rather than intra-halo departures from 
a smooth functional form (i.e. "bumps"). For small radii 
(r ^ r-2), the Gaussian PDF of halo spherical shell densities 
is well-matched by assuming that each halo is described by a 
deterministic Einasto radial density profile whose concentra- 
tion is drawn from a Gaussian distribution (see Fig. This 
implies that the distribution of halo concentrations fully de- 
termine the PDF of the smooth density component, and 
provides additional support that the form of the concentra- 
tion PDF is Gaussian. More specifically, the effect of radial 
density "bumps" within the halo cannot be larger than the 
effect of measurement uncertainty of individual halo concen- 
trations. A log-normal concentration distribution is unable 
to reproduce the density distribution shape or peak. Unsur- 
prisingly, an assumption of an NFW profile does not match 
the density peak, although it is able to produce the Gaussian 
shape and width. 

At larger radii (r XL f-2), the distribution remains 
Gaussian but is wider than implied by the distribution of 
concentrations, possibly due to increased scatter introduced 
from substructure or from a lower degree of relaxation in 
the outskirts of these recently formed cluster-size halos. The 
distribution of densities in our data is described to (~ 10%) 
accuracy by the following function: 

logiodp = 1.144 login (p) - 1-389, (6) 

where (p) and (jp are found via a Gaussian fit to the densi- 
ties measured from the halo sample at a given radius. The 
agreement between the fit and the simulation data, shown 
in Fig. [S] over more than three orders of magnitude in den- 
sity is interesting. However, we do not advocate that this 
is a universal function; further work is required to deter- 
mine whether this relation between density and its scatter 
remains valid for lower halo masses and different redshifts. 
The moderate flattening of Op at low density reflects the 
relative widening of the density PDF at large halo radii. 

It is instructive to consider also the distribution of halo 
density profile slopes. In Fig. |6l we show the distribution of 
the logarithmic radial slope of the density at various radii. 
The PDF of din p/dlnr is well-described by a Gaussian dis- 
tribution, except perhaps at large radii where substructure 
or other effects appear to result in wider than Gaussian 
tails. The mean of this distribution at each radius is consis- 
tent with the Einasto profile with a — 0.19. Note that the 
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Figure 2. Left panel: Density profiles for 100 random halos from the sample. Points with error bars denote the mean and the Icr scatter 
of the distribution for the full halo sample. Right panel: Logarithmic density slope for the 100 largest halos. Slopes are normalized to 
the scale radius to remove the effect of scatter in halo concentrations. An Einasto profile with a = 0.19 is a good match to the mean 
simulation slope. 




Figure 3. Halo-to-halo distribution in density at various radii 
relative to r200 in the full halo sample. 

Einasto profile with fixed a implies zero scatter in the PDF 
of halo slopes. The radii in the right panel of Fig. |6] are cho- 
sen such that they should have identical logarithmic slopes, 
assuming the Einasto profile, for the mean halo concentra- 
tion of 4.55. For most radii, the width of the distribution of 
slopes are remarkably similar whether measured in terms of 



Figure 4. Measured density distribution versus density distri- 
bution predicted by a deterministic Einasto profile and a Gaus- 
sian concentration distribution of width crc200 = 1.28, centered 
on C200 = 4.55, to match the measured values of our sample. 
Log-normal concentration prediction assumes the best-fit to our 
sample of o-jogjgj. = 0.121 with median c = 4.66. Gaussian con- 
centrations match the shape and peak position of the density 
distribution better than log-normal concentrations. 
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Figure 5. Gaussian Op as a function of mean density at fixed 
"r/fvir (black crosses) and fixed r/r200 (blue squares) for our halo 
sample. The density distribution, P{p), is well-described by a 
Gaussian distribution at all radii. Line is a fit given by Eq. [B]!. 



r_-2 or in terms of r2oo, apart from the innermost plotted 
radius. This suggests that intra-halo "bumps" rather than 
halo concentration is the major contributor to scatter in the 
slope PDF. In fact, Eq. [T] implies that the distribution in 
density slopes at fixed r/r2oo due to the concentration dis- 
tribution should be ^0.1, which is much smaller than the 
actual spread which ranges from O.IT <C f^dinp/dinr 

< 0.34 

at the fixed r/r2oo values shown. At the smallest radius of 
r_2/10, the slope distribution width is likely dominated by 
errors in the concentration measurement (left panel). Pre- 
sumably for this reason, the PDF of the slopes is narrower 
when considered with respect to r2oo at the smallest radii. 
At each radius, we show an estimate of the uncertainty in the 
slope measurement, based on poisson noise from the average 
number of particles per radial bin. The slope uncertainty 
is significantly smaller than the measured PDF, except at 
the smallest plotted radius. Because numerical problems are 
most difficult to overcome at small radii, poisson uncertainty 
could underestimate the true error. We thus cannot rule out 
the possibility that the broadened PDF at small radii may 
have numerical origins. Indeed, increasing the minimum halo 
mass of the sample by a factor of 5 narrows the slope PDF at 
the smallest radii such that this effect is significantly smaller. 



5 EFFECTS OF PROFILE SCATTER ON 
DARK MATTER ANNIHILATION 

In this section we discuss the effect of the profile scatter on 
the expected signal in 7-rays (or other byproduct) from dark 
matter annihilation in halos. 

In general, the total gamma-ray luminosity from a halo 



of mass M is given by the volume integral of the square of 
its mass distribution as 

C^^^j\Ur)d\ (7) 

where M-^ is the mass of the dark matter particle, and (tju) 
is the thermal average of the annihilation cross section. Par- 
ticle physics enters through the mass of the dark matter 
particle, and through its total annihilation cross section. As 
an example, in order to demonstrate the effects of density 
distributions on the annihilation flux, we consider a dark 
matter particle with mass = 400GeV, with a total anni- 
hilation rate to bh quarks given by (cr«) — 10~^^cva? /s. We 
assume these values throughout the rest of this manuscript, 
and note that in general, the assumed dark matter particle 
properties affect the normalization of the results presented, 
and not the shape of the distributions. 

The "dark luminosity" of a halo is determined by the 
distribution of its mass (see Eq. [7|. We assume that the 
mean distribution of dark matter in halos is described by 
the Einasto profile, from Eq[l}(4l with a — 0.1645 as given 
in Gao et al. (2008. ) (and we ignore any mass and redshift 
dependence of a). For the remainder of this section, we focus 
on the dark luminosity of the smooth component in a halo 
and thus we ignore the presence of substructure and any 
associated "boost" that may contribute to the annihilation 
luminosity as defined in Eq. [7] 

5.1 The flux distribution as a function of mass 

We first confirm that we can capture the PDF of the lo- 
cal dark matter annihilation volume emissivity (in spheri- 
cal shells). This is the quantity that will be integrated to 
compute total halo annihilation luminosity. In Fig. [T] we 
compare the measured distribution of normalized differen- 
tial annihilation luminosity per logarithmic radial interval 
(47rp^(r)r^), shown here at r = 0.032r2oo as an example, 
with that from an Einasto profile with the mean concentra- 
tion and Gaussian scatter of the halo sample (Eq. [71 . We 
plot this distribution in units of mean density and r2oo so 
that the quantity is independent of halo mass. The excel- 
lent agreement of measured and predicted differential anni- 
hilation luminosity implies that the concentration distribu- 
tion with the assumption of an Einasto profile is sufficient 
to estimate localized (in radius) dark matter annihilation 
luminosity. However, we have yet to establish that the lo- 
calized annihilation strength can be integrated to yield the 
correct halo annihilation luminosity. Correlations of density 
with radius could result in large scatter in annihilation lumi- 
nosity from halo to halo. For example, halos that happen to 
have enhanced density over some extended range below r_2, 
where the annihilation rate is larger (albeit within a smaller 
volume), could have enhanced annihilation luminosity ver- 
sus halos with smoother profiles and similar concentrations. 
Although the results of §|4]iinply that intra-halo radial corre- 
lations are less important than the scatter in concentrations 
on local density, the p'^-dependence of annihilation could 
lead to a larger impact on halo annihilation rates. 

In order to assess the origin of the scatter in halo anni- 
hilation luminosities and the contribution due to the scatter 
in concentrations, we compute the dark luminosity of each 
simulation halo, and compare with the expected luminosity 
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Figure 6. Distribution in halo slopes at various radii (solid lines) relative to the scale radius (r_2 = ?'20o/c20o) (left), and relative to 
f'200 (right). Dotted curves are best-fit Gaussians. Vertical lines denote the slope of the Einasto profile with a = 0.19. In the right panel, 
the Einasto profile slopes are shown for the mean halo concentration of 4.55. Slopes are reasonably fit by a Gaussian with an average 
consistent with the Einasto profile and a width that is narrowest near . Horizontal lines indicate an estimate of the slope measurement 
uncertainty. 
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M[h-iMg] M[h-iMg] 

Figure 8. The inferred dark matter luminosities in units of 10™GeV^cm~'^. Each dot represents a dark matter halo, while the vertical 
bars represent the 68 percentile of the distribution in each virial mass bin, whose width is denoted by the horizontal bars. Left panel: 
Luminosity computed directly from the density profile extracted from the simulation for each halo. Right panel: Expected luminosity 
computed from fitting a concentration to each halo and assuming an Einasto density profile. 
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Figure 7. Normalized dark matter annihilation "luminosity" 
per logarithmic radial interval (47rp^ (r)r'^) for the full halo sam- 
ple. The prediction assumes an Einasto profile with a Gaussian 
concentration scatter, as in Fig. |4l 



5.2 Cosmological 7-ray background 

We now consider the contribution of the scatter in densities 
to the cosmological 7-ray background: the annihilation flux 
integrated over all halos at all masses and redshifts. We are 
interested in the effects of the non-universality of profiles to 
the expected gamma-ray background. 

We compute the gamma-ray background as 



dEdAdtdO. 



1 c 





dn{z) 



(av) 



dN^ 



h(z) dMdV 



Ml d[E{l + z)] 
pliir)d^rdM dz, 



(8) 



where c is the speed of hght, Ho — 0.7 i s the present value 
of the Hubble constant, and h{z) = y^O,M{i + z)'-^ + Qa- 
We assume that Qm = 0.27 and SIa ~ 0.73. The quantity 
dN^ /d[E{l + z)] is the spectrum of the emitted photons at a 
source energy of E{1 + z), and the mass function of objects 
of mass M at redshift z is dn (z) /dMdV . We use the mass 
function of iReed et aTl (|2007| ). 

We assume that the annihilation proceeds to a bb 
quark final state, and that the distribution in the num- 
ber of gamma-rays emitted per source energy inter- 
val Ea are described by the fu nctional form given in 
iBergstrom. Ullio fc Bucklev I l|l998l ). namely, 



given its measured concentration. We show in Fig. [8] the dis- 
tribution of halo dark luminosities computed directly from 
the simulation halo density profile (left panel) together with 
the prediction of the same quantity from an Einasto profile 
using the individually measured concentration of each halo 
(right panel). We bin the data in mass, and determine the 
mean and 68 percentile of the distribution. The measured 
and predicted luminosities and 1 a scatter agree well; this 
suggests that radial correlations in density should not pre- 
vent accurate estimation of the annihilation luminosity of a 
halo. 

Thus, the origin of the distribution of luminosities at 
each mass bin is the distribution in concentrations, which 
correlate with formation time, albeit with large scatter (see 
e.g. iNeto et aT~ll2007l ). The correlation between concentra- 
tion and annihilation luminosity can be modelled in the fol- 
lowing manner for the smooth density component of dark 
matter halos. The normalization of the profile is propor- 
tional to p-2 ~ CjQo, where 7 ~ 3 for C200 ^ 1, and 
r_2 ^ c^oo- follows that, roughly speaking, the luminos- 
ity scales as L ~ pL2'''-2 ~ '^200^' leading to L ~ cloo for 
C200 S> 1, although for concentrations typical of our clus- 
ters L ~ C200 ■ Note that it is difficult to distinguish between 
the subtle differences between a Gaussian distribution of 
concentrations and a log-normal one from Fig. [S] In addi- 
tion, we note that the mean of the luminosity distribution 
at each mass bin is roughly proportional to the mass of the 
halo. This is to be expected as the luminosity of a dark mat- 
ter distribution that is described by a two-parameter profile 
(e.g., NFW, and/or Einasto) is L ~ p?.2'"-2 ~ , where 
/3 ~ 1 because the dependence of concentration on mass is 
relatively weak. 



dErM-Xw^) '"'"''Y wA' 

We estimate the cosmological gamma-ray background 
for three different halo concentration distributions. First, we 
assume an Einasto functional form of the density as a func- 
tion of radius and a one-to-one d ependence of c o ncent rations 
on mass and redshift as given in lMaccio et al.l(|2007l '). How- 
ever, halo concentrations exhibit a distribution at fixed halo 
mass and redshift (see ii l4.1[) . As such, we consider a second 
case in which halo concentrations instead follow a Gaussian 
dis tribution, given by a g = 0.283c2oo, where C200 remains 
the'Macci o et al."1 (|2007l ) concentration as a function of mass 
and redshift, and the Einasto profile still describes the den- 
sity distribution. This value for the Gaussian width corre- 
sponds to the fractional width found in our halo sample. 
Finally, we consider a case where the concentration distri- 
bution is log-normal about the mean concentration value, 
with dispersion given by aLN = 0.121, and all other aspects 
of the background calculation remain the same. 

In Fig. [9l we show the expected cosmological gamma- 
ray background for the three different distributions of dark 
matter. As expected, in the presence of a spread in the distri- 
bution of concentrations, the annihilation flux is increased 
relative to the case where there is a one-to-one mapping 
between concentration and mass. We find that a log-normal 
distribution would give rise to approximately a 10% increase 
in the annihilation flux relative to our preferred Gaussian 
distribution of concentrations (which is only a few percent 
higher than the simple case of no distribution in concentra- 
tions). 
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Figure 9. The effects of a Gaussian {short- dashed) and log- 
normal (long-dashed) distribution of concentration parameters 
to the cosmological extragalactic dark matter annihilation back- 
ground. The solid curve represents the background computed by 
assigning a fixed c(M) relationship without a spread. In all cases, 
the background is assumed to be due to a 400 GeV WIMP an- 
nihilating with a cross section of 10~'^^cm^s~^ into a bb quark 
pair. 



Figure 10. The probability distribution function of gamma- 
ray flux along a line of sight at 30° {short- dashed), 90° {long- 
dashed) and 180° {dot-dashed) degrees from the Galactic center 
for a Gaussian concentration distribution {thick curves) and a log- 
normal concentration distribution {thin curves). Note the smaller 
width of the distribution at high Galactic angles relative to the 
inner parts of the halo. Particle physics parameters are same as 
Fig.H 



5.3 The annihilation flux due to the smooth 

distribution of dark matter in the Milky Way 

We now discuss the impact of the halo density probabil- 
ity distribution function on the Milky Way annihilation flux 
along different lines of sight. The expected flux at a partic- 
ular angle ijj with respect to the Galactic center depends on 
the distribution of dark matter densities along that line of 
sight. As that is an outcome of the particular concentration 
of the Galactic Halo, drawn from a distribution of possible 
concentrations, there is a distribution of expected fluxes for 
each line of sight. 

The line of sight flux at an angle ^j) with respect to the 
Galactic center can be written as 



dEdAdtdO, 



1 {av) dN-y 
4^ Ml ~dE 



where €max = [cos ijj -\- a/ {Rmw /dQ^ 



' tp\ and 



yJd'^-\- P — 2do^ cos ■(/). We take the distance of the Sun 
from the Galactic cen ter to be = 8.5kpc (consistent with 
iGillessen et al. Il2009( ) , and the radius of the Milky Way halo 
R200,MW ~ 250kpc, which implies M200 = 1.8 x IO^^Mq ari d 
M^ir ~= 2.xlO^^M0 f consistent withe.g. lCuo et al. l20l"oh . 

In Fig. 1101 we show the expected flux distribution at 
various angles with respect to the Galactic center computed 
using Eg. 1101 We calculate the flux distribution for two cases, 
first where the distribution of concentrations is Gaussian, 
and second where the distribution of concentrations follows 
a log-normal distribution. We find that for a log-normal dis- 
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Figure 11. The dependence of the peak of the distributions 
(Gaussian/log-normal) to the angular distance of the line of sight 
with the Galactic center. Log-normal distribution in concentra- 
tions leads to a slightly higher peak in the flux probability distri- 
bution function. 
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Figure 12. The dependence of the width of the distributions 
(Gaussian/log-normal) to the angular distance of the line of sight 
with the Galactic center. Log-normal distribution in concentra- 
tions leads to a slightly narrower flux probability distribution 
function. 



tribution of concentrations the width of the flux distribu- 
tion along a line of sight is slightly narrower, while at the 
same time, the mean of the distribution is slightly higher. 
This is to be expected as the annihilation rate is sensitive 
to concentration parameter. The high concentration tail of 
the log-normal concentration distribution contributes to its 
higher mean flux, while the more extended low concentra- 
tion range from the Gaussian distribution manifests itself 
into a broader distribution of fluxes at each particular an- 
gular Galacto-centric distance. 

It should also be emphasized that the shown distribu- 
tion functions are uncorrelated, while in reality, because a 
single concentration must be defined for the Halo, there are 
correlations between the flux at adjacent angular bins, and 
anti-correlations between angular values close to zero and 
180 degrees. A higher concentration halo has relatively more 
mass near the center and less near the outer parts. Thus, 
highly concentrated halos will have higher fluxes with re- 
spect to the distribution function toward the Galactic center, 
and will have relatively smaller fluxes toward the Galactic 
anti-center. 

We now quantify the expected angular dependence of 
the peak and width of the flux probability distribution func- 
tion. The peak of the flux distribution is expected to be 
smaller at high angular distances from the Galactic center. 
This is a natural consequence of the centrally concentrated 
spherical distribution of dark matter in a halo, and enables 
a measure of the underlying density profile of the halo. In 
Fig. 111! we show the angular dependence of the peak fiux. 
A good fit (within few percent) to this function is obtained 
by a double power-law as: 



Gaussian 


-13.70 


-0.39 


-1.95 


3.27 


-0.12 


0.98 


Log — normal 


-13.48 


-0.40 


-1.96 


3.22 


-0.11 


0.87 



Table 1. Fitting parameters of Eg. 1111 fc 1121 for the peak and 
width of the gamma-ray annihilation flux distribution function 
along a Galactic line of sight. 



where the parameters ai, a2, 03 and 04 are given in Table [T] 
The width of the distribution is smaller at large angles from 
the Galactic center. This is to be expected because the ef- 
fects of concentration are more apparent in the inner regions 
of the halo. At radial distances r 2> r-2, the changes in 
the dark matter density due to different halo concentration 
values are smaller and therefore the flux distribution is nar- 
rower. In Fig. 1121 we show the expected angular dependence 
of the width of the flux distribution. We find that a function 
of the form 



(12) 



In <1> = ai 4- 02 In ?/) -)- 03 ln(a4 + V") 



(11) 



is a good (within few percent) fit to the angular dependence 
of the width of the distribution function. The quantities as 
and ag are given in Table[T] It is important to note here that 
for large radii (r ^ f-2) the halo density PDF that we mea- 
sured in § 14.21 is larger than inferred from the concentration 
scatter, which implies that the values of flux uncertainty 
at large Galacto-centric angles may be larger than our es- 
timates. However, because the Solar Radius is well within 
r_2, the flux should be dominated by the density distri- 
bution within r_2, even toward the Galactic anti-center, so 
effects of large radii scatter on the flux PDF should be small. 



5.4 Other applications 

5.4-1 Cosmological Uncertainty of the Local Dark Matter 
Density 

Applying our results, as before, to a Milky Way mass of 
M200 ~ 1-8 X IO^^Mq, a concentratio n of C200 = 5.97 (us - 
ing the concentration-mass relation of lMaccio et al. ||2008| ). 
and a Solar radius of 8.5 kpc implies a Solar radius to- 
tal matter density of 0.262t[;j^GeV cm if the concen- 
tration PDF is assumed to be Gaussian with width pro- 
portional to our values of CTc/c = 0.283. Assumptions 
of a log-normal distribution of width gln = 0.121 re- 
sults in only minor changes for a local density range of 
0.262^0343 GeVcm"^, which can be compared with sev- 
eral o bservational estimates. iBergstrom. Ullio fc Bucklev I 
(jl998l ) find an allowed ra nge of [0.2-0.8] GeVcm "^, and 
a more recent work by iWeber fc de Boerl l|2010t ) find an 
acceptable range of [0.2 -0.4] GeVcm"'^. However, t ighter 
constraints are found bvlW^idrow. Pvm fc Dubinski I (|2008D 
and ICatena fc Ullio I (|2009l ). who utilize a varietv of dy- 
namical observables to estimate, respectively, 0.304 ± 0.053 
GeVcm"^, and 0.385 ±0.027 GeVcm"^ for the dark matter 
density. These estimates are somewhat larger than our cos- 
mological range, which hints at the possibility that the Solar 
radius dark matter density has been enhanced by "adiabatic 
contraction" in response to baryon cooling. 
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5.4-2 Implications for Halo Stacking 

Our results have implications for many astrophysical ap- 
plications that depend upon the mass distribution within 
halos. One such example is the technique of "stack- 
ing" halos to improve signal to noise, commonly employed 
in simulations and observations. In one such applica- 
tion, large numbers of simulated halo density profiles 
are stacked to measure the mean density profile to hig h 
precision (e.g. ICao et al. I l2008l : iHavashi fc White I l200a ). 
From an observational perspective, stacking many halos 
of similar mass greatly reduces the noise in, for exam- 
ple, weak lensing de terminations of halo mass profiles and 
concentrations fe.g. [ Mandelbaum. Seliak. fc Hirata I l2008l : 
iMandelbaum et al IbOoi T 

Our results support the viability of halo "stacking". 
Due to the fact that the distribution of densities at fixed 
radius is Gaussian and the distribution of concentrations is 
also Gaussian, a stacked density profile is indeed an accu- 
rate representation of the median profile. This is particularly 
convenient in that it allows the mean halo profile to be pa- 
rameterized into an analytic form without bias, and allows 
unbiased stacking of observational mass profiles. In Fig. 1131 
we have verified that the cumulative cluster mass distribu- 
tion also remains unbiased. We compare the Einasto form 
of a 3-dimensional enclosed mass profile of C200 = 4.55 with 
the same quantity for a mock stacked halo drawn from a 
Gaussian distribution of concentrations of mean C200 ~ 4.55 
and scatter ac = 0.283c2oo, and find agreement to better 
than 1%. This implies that mass profiles determined by lens- 
ing studies should be free from biases associated with halo 
stacking. If, instead, the density distribution had been log- 
normal, then a stacked halo would have been biased toward 
higher cumulative masses at small radii, although not by 
more than a few percent. 




0.01 0.1 1 

Figure 13. Enclosed mass in spherical shells for mock stacked 
halo samples with Einasto profile at C200 = 4.55 plus Gaussian 
concentration distribution and Log-normal concentration distri- 
bution. C200 = 3.27 and C200 = 5.83 denote the 1 — a range in 
concentration distribution. Bottom panel shows the ratio of each 
mass profile with respect to mean concentration, C200 = 4.55, 
Einasto profile. Level of agreement with mean concentration in- 
dicates robustness of M < r in stacked halo samples against bias 
caused by halo-to-halo scatter in concentrations. 



6 LIMITATIONS OF THIS WORK 

In this work, we do not consider other potential impor- 
tant factors to the annihilation rate, which could possi- 
bly be dominant over the halo-to-halo scatter associated 
with hierarchical structure formation (via the distribution 
of concentrations) that we have examined. Among them is 
the amount of substructure in small subhalos and streams 
(the "clumpy" dark matter component), whose contribu- 
tion can boost the annihilation rate relative to that of a 
smooth halo, and will ad d to the uncertainty in the lo- 
cal dark matter density (e.g .'Kamionk owski fc Koushiappas I 
l2008l : iKamionkowski et ahfeoiO. . Vogelsberger et al. \\20oW . 
We also ignore any gravitational coupling that the dif- 
ferential evolution of baryonic halo component may have 
on the dark matter halo structure. Baryon infiuences may 
include gas cooling; this could cause the dark matter 
halo to res pond to the deeper pot ential by "adiabatic 
contraction" (iBlumenthal et ahl Il986l) . Although these ef- 
fects may be important, they are beyond the scope of this 
study. 

Our application toward annihilation rates in the Galaxy 
relies upon the assumption that the behavior of the halo- 
to-halo profile scatter is similar for galaxies and clusters, 
namely that the probability distribution of the mean den- 



sity in radial shells is always described by the Einasto pro- 
file with a Gaussian distribution in halo concentrations. 
The assertion that the distribution of halo concentration 
remains universally Gaussian, while speculative, is sup- 
ported by the fact that the logarithmic width and shape 
of the concentration dist ribution has weak or no mass or 
redsh i ft dependence (e .g. Bullock et a . _ 200 ll : iNeto et al. I 



2OO7I: iGao et al. I l2008l: iMaccio et al. Il2007l l. Additionally, 
Knollmann. Power fc Knebe I (|2008 ll used scale-free simula- 



tions to show that the scatter in halo profile concentration 
and density slope has little dependence on matter power 
spectral index (which varies with halo mass) over a range 
bracketing well beyond the effective spectral indices of clus- 
ters and galaxies. However, for small halo masses, the range 
in formation redshifts is much broader than for clusters, 
which have formed recently or are still undergoing the virial- 
ization and relaxation processes associated with formation. 
Given the link between halo concentration and formation 
redshift, this could lead to a more significant high concen- 
tration tail for lower mass halos, and it could be argued that 
this should favor a log-normal distribution in halo concentra- 
tions, rather than the Gaussian distribution for our clusters. 
Future studies are warranted utilizing a wider range in halo 
masses. 
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7 CONCLUSIONS 

The probability distribution function of dark matter within 
halos, as we have explored in this work, provides some basis 
for interpreting both indirect and direct dark matter de- 
tection experiments in a cosmological context. Constraints 
upon the dark matter density, particle mass, or the self- 
annihilation cross section depend on the probability distri- 
bution function of dark matter. 

Our results indicate that halo concentration is the pri- 
mary cosmological contributor to the dark matter PDF. 
This implies a particular correlation between the local dark 
matter density, relevant for direct detection efforts, and the 
dark matter density in the direction of the Galactic cen- 
ter (and elsewhere), applicable to indirect detection exper- 
iments. The effect of halo concentration should thus be a 
crucial factor in verifying the consistency of dark matter 
density constraints made from multiple dark matter detec- 
tion techniques. Ultimately, dark matter signals might be 
able to test the validity of the ACDM cosmological model 
through estimates of the dark density at differing locations 
within the Milky Way halo, and perhaps also within other 
halos. 
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